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Abstract 

We establish a new algorithm that generates a new solution to the Einstein 
field equations, with an anisotropic matter distribution, from a seed isotropic 
solution. The new solution is expressed in terms of integrals of an isotropic 
gravitational potential; and the integration can be completed exactly for par- 
ticular isotropic seed metrics. A good feature of our approach is that the 
anisotropic solutions necessarily have an isotropic limit. We find two exam- 
ples of anisotropic solutions which generalise the isothermal sphere and the 
Schwarzschild interior sphere. Both examples are expressed in closed form 
involving elementary functions only. 
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1 Introduction 



Numerous models of static perfect fluid spheres, in the context of general relativity, 
have been constructed in the past because these are first approximations in building 
a realistic model for a star. Lists of exact solutions to the Einstein field equations 

nnn 

modelling relativistic perfect fluid spheres are given in several treatments I2L L| • 
In these works it is assumed that the matter distribution is isotropic so that the 
radial pressure is the same as the transverse pressure. A strong case can be made 
to study matter distributions which are anisotropic in which the radial component 
of the pressure is not the same as the transverse component. Anisotropy should not 
be neglected when analysing the critical mass and the redshift of highly compact 
bodies, and it is important in modelling boson stars and strange stars. Consequently 
anisotropy has been studied extensively in recent years by a number of researchers 



within the framework of general relativity 



recent years by a 



hi 



Most solutions of the Einstein field equations with isotropic matter have been 
obtained in an ad hoc approach. Recent investigations have attempted to generate 
isotropic models using a systematic and algorithmic approach Q, Q|. How- 



ever there exist very few analo gou s results for generating anisotropic models. In 
this regard Maharaj and Chaisi ly] have established an algorithm that produces a 
new anisotropic solution from a given seed isotropic line element. Here we present 
a different algorithm that generates a new anisotropic solution. The present algo- 
rithm has a simpler form and involves fewer integrations, and it is consequently 
easier to apply. A desirable physical feature of our approach is that these solutions 
have an isotropic limit; we expect that gravity acts to eventually isotropize matter 
in the absence of other external forces. Note that many of the exact solutions found 
previously remain anisotropic with no possibility of regaining an isotropic matter 
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distribution in a suitable limit 
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The main objective of this paper is to demonstrate that it is possible to generate 
anisotropic solutions from a given isotropic solution. To achieve this we need to 
complete an integration which we show is possible for particular metrics. In §2 we 
provide the fundamental field equations for isotropic and anisotropic matter. The 
field equations are presented as a first order system of differential equations. The 
algorithm that produces a new anisotropic solution is described in §3. As a first 
example we use the isothermal model to produce a new anisotropic solution in §4. 
As a second example we use the Schwarzschild interior model to produce a new 
anisotropic solution in §5. In both examples the solution can be given explicitly in 
terms of elementary functions which makes it possible to study the physical features. 
We briefly study the behaviour of the anisotropy factor in both examples. 



2 Field equations 

We utilise a form of the Einstein field equations in which only first order derivatives 
appear. This representation assists in simplifying the integration process as pointed 
out by Chaisi and Maharaj 19j| whose notation and conventions we follow. The line 
element for static spherically symmetric spacetimes is given by 

ds 2 = -e u dt 2 + e x dr 2 + r 2 (d6 2 + sin 2 6d<j) 2 ) (1) 

where v(r) and A(r) are arbitrary functions. The energy-momentum tensor, for 
nonradiating matter, for isotropic distributions has the form 

T ab = {fi + p)u a u b + pg ab (2) 



3 



where /i is the energy density and p is the isotropic pressure. These are measured 
relative to the comoving four-velocity u a = 6~ u ^ 2 5q. We define the mass function as 



m ( r ) = o / x 2 /i(x)dx (3) 

Consequently M = m(R) is the total mass of a sphere of radius R. The Einstein 
field equations are equivalent to the system 

e~ A = 1 (4) 

r 

r(r — 2m)v' = pr 3 + 2m (5) 
(^ + p)u' + 2p' = (6) 

where we have used 

The energy-momentum tensor for anisotropic matter which is not radiating has 
the form 

T ab = (/j + p)u a u b + pg ab + n ab (7) 

The quantity 7r ab = \/3S(r) (c a c b — lh ab ) is the anisotropic stress tensor; the space- 
like vector c a = e~ A//2 5" is orthogonal to the fluid four-velocity u a = c~ u / 2 5q and 
\S(r)\ is the magnitude of the stress tensor. The Einstein field equations, with the 
metric (p^ and the matter content ( | 7 j) . can be written in the form 

\ 2m . . 

e- x = 1 8 

r 

r(r — 2m)v = p r r 3 + 2m (9) 
4 

{/I + Pr ) V' + 2p' r = (Pr-P±) (10) 

r 

for anisotropic matter distributions. The radial pressure p r is distinct from the 
tangential pressure p±. It is convenient to write p r and p± in the form 

p r = p + 2S/ V3, p± = p - S/y/3 

where S provides a measure of anisotropy. Note that for isotropic matter p r = p± = p 
and we regain (jUj-dHJ). 



3 The Algorithm 



In this section we establish a procedure for generating a new anisotropic solution 
of the Einstein field equations from a specified isotropic solution. We start by 
considering the Einstein field equations (0}-© with isotropic matter distribution. 
We assume that an explicit solution to (0)-© is known where 

(u,X,m,p) = (v , \ ,m ,p ) (11) 

and functions z/ , Ao, m and po are explicitly given. Then the equations in (JH)-(JB1) 
are satisfied and we can write 

= l-^o (12) 
r 

r(r — 2m )u' = p r 3 + 2m (13) 
(^T + P°) "0 + 2 Po = (14) 

We next consider the Einstein field equations (fKjl- lfTUJ) with anisotropic matter dis- 
tribution and seek an explicit solution. To this end we propose the possible solution 

( a(r)\ 
(v,\,m,p r ,p ± ) = \v ,\ + x(r),m + y(r),p + a(r),p (15) 

where (z/ , A , m ,p ) are given by (fTTj) and x,y and a are arbitrary functions, and 
we have set a = 2S/V3 for convenience. Then the system (P)- (fTT)|) becomes 

-(Ao-Hr) i 2m o + 2 V 



e -(Ao+x) = 1 (i 6 ) 

r 

r(r — 2m — 2y)u' Q = p r 3 + ar 3 + 2m Q + 2y (17) 

2m' + 2y' \ . n , n . 6 

+ po + a J i/ + 2p' + 2a' = — a 
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The systems and lead to 



./■ - In { 1 - — e A ° } (19) 



ar 3 



2(l + ri/ ) 



(20) 



2y' \ . . 6a . . 

■£ + a)i/ + 2a' = -- (21) 

We need to integrate ()21)1 to find the function a. The remaining functions x and 
are defined in terms of a. Two cases arise: a = and a^O. If a = then (|19 |1 -(|21 )1 
has the solution 

(a,x,y) = (0,0,0) (22) 

The trivial solution (j22j) corresponds to the isotropic case. Thus this algorithm 
regains the isotropic solution in the appropriate limit. If a ^ then we can eliminate 
y from (|2~Tj) to get 

2i/ n ( a'r 3 3ar 2 ar 3 v' n + rv'l 1 , , 6a 



r 2 [ 2(1 + ri/ ) 2(1 + ri/ ) 2 (1 + ru' ) 2 J u r 
This differential equation can be written as 

a' v' n f„ ^o + ^qI ^ ' ■ 6\ f\ + rv' ( 



< + Z hrhdM (23) 



a 2 + ru' Q \ I + ru' Q J \ r J \2 + rv' Q 

after some simplification. On integration (}23|) leads to 



In a = J Q + lnfc (24) 
where In /c is a constant of integration, k ^ 0, and we have set 



2 + rv' Q \ l + ru' J \° r J \2 + ru' 
We can write ()24|) in the compact form 



a = ke Ja (25) 
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Equations (JUJ), (J27?j) and correspond to anisotropic matter. 



Thus if given a known isotropic solution (jllj) we can generate a new anisotropic 
solution (fT5|) where 



a = ke Ja (26) 
x = -ln(l-^e A °) (27) 



r 

„3 



ar 

V = -T777-T^77T C« 



2(1 + ri/ ) 



and the integral J a is given by 



The integration in J a can be explicitly performed as i/q is specified in the isotropic 
solution (JTTJ). Note that (|2l)j) - (j2*8j) applies to both cases a = and q^O. If a = 
we can set A; = and regain the isotropic result (|22jl . When a ^ then 7^ and 
we regain the anisotropic equations (JTHJ), (jSOj) and (J2*H|) . 

It is remarkable that our simple ansatz leads to a new anisotropic solution of the 
Einstein field equations. This is subject to completing the integration in J a ; clearly 
this is possible for particular choices of the isotropic function i/q. We demonstrate two 
examples of anisotropic solutions for familiar choices of z/ in the next two sections. 
The algorithm that we have generated in this paper is easy to apply as there is only 
a single integration to be performed unlike the earlier algorithm of Maharaj and 
Chaisi 3] which is more complicated and involves further integrations. We believe 
that new anisotropic solutions that arise from our procedure are likely to produce 
realistic anisotropic stellar models. We emphasise that a desirable feature of our 
approach is that our models contain an isotropic limit which is often not the case 
in other approaches. 
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4 Example 1 



As a first example we demonstrate the applicability of the algorithm in §3 by gener- 



ating anisotropic isothermal spheres. The line element for the isothermal model 
has the form 



21] 



ds 2 = -r^dt 2 + (l + -jY^y) dr2 + r2 ( d ° 2 + sin2 6d< f) 



(29) 



where c is a constant. The corresponding isotropic functions for ()29j) are given by 

4c , , f 4c 1 2cr 4c 2 / 1 
In r, In < 1 H > , ; — , ; — - — 

1 + c ' I 1 + c J '4c+ 1 + c 2 ' 4c+ 1 + c 2 \r 2 



^0 = J" , , „N9 ( 3 ) ( 31 ) 



The energy density function has the form 

4c + (1 - r)- V/'- 
Hence ()30|) and (}3T|) imply that 

Po = cfio (32) 

which is a linear barotropic equation of state. Isothermal spheres with the density 
profile n oc r~ 2 and the equation of stateJ22j) appear in a variety of models for 
both Newtonian and relativistic stars jljj, They have been studied extensively 
in astrophysics as an equilibrium approximation to more complicated systems which 
are close to a dynamically relaxed state j^. 

With the isotropic functions (jHOJ) we can evaluate the integral J a in (j2Hj) and we 
find 

, ,'3 + 14c+ 19c 2 \ dr 3 + 14c + 19c 2 , , , 

J n = — — — = — In r + In k 

' 1 l + 4c + 3c 2 J r l + 4c + 3c 2 
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which leads to the expressions 



_ 3+14c+19cf 

a = kr i+4c+3c^ 

. J\ 4 C +(1+C 2 ) _ l+6c+13c; 

x = - In < 1 + k- 77 —r i+ic+sc 2 

\ (l + c)(l + 5c) 

k ( 1 + c 

2/ = 



2 \l + 5c / 

Consequently we obtain the new line element in the form 



4c 



4c \ ( , 4c + (1 + cr _i+6c+i3c2 

ds 2 = -r— dt 2 + 1 + - _) M+fc- y- ^— r i+4c+3^ dr 2 

(l + c)V V (l + c)(l + 5c) 

+r 2 (dfl 2 + sin 2 Odcj) 2 ) (33) 



-i 



and the matter variables have the analytic representation 

2cr k ( 1 + c \ 2c+io e 2 

m = ; r- r i+ 4c +3c 2 (34) 

4c + (1 + c) 2 2 Vl + 5cy V ; 

1 4C 2 , 3+14 c +19c 2 

4c + (1 + cr 

1 4C 2 3+14c+19e 2 

rj = r i+4 C +3c^ (35) 

P± r 2 4c+(l + c) 2 2 1 ; 

The isotropic isothermal sphere model ()29|) produces the anisotropic isothermal 
sphere model (|3*3*|) utilizing our algorithm. With the parameter values k = 0, we 
regain the conventional isothermal sphere. 



The degree of anisotropy is 

k r- 3+14c+19c 2 

S = -V3r i+4 C +3c 2 (37) 

Mathematical^ was used to graph the anisotropy factor (pT7|) . The plots are as 
shown in Figures ^ and El for particular values of the parameters shown. The 
anisotropy factor S is plotted against the radial distance on the interval < r < 1. 
There is a singularity at r = that has been carried over from the other dynamical 
and metric functions. However because constants k and c can be picked arbitrarily, 
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the pair can be chosen such that S(r) is monotonically decreasing or increasing. The 
physical considerations of a problem may lead to the choice of one profile as opposed 
to the other; for example the S(r) profile in Figure Q may be preferable where a 
stellar body with vanishing anisotropy as one moves from the center of the body to 
the boundary is considered. The S(r) profile in Figure El could be chosen over the 
one in Figure Q for boson star models as proposed by Dev and Gleiserj^]. The fairly 
simple behaviour of the S(r) in these plots shows that a more extensive physical 
analysis of the solutions is possible, which will be carried out in future work. 
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Figure 1: Decreasing S(r) for anisotropic isothermal sphere; c = 1, k = .01 




Figure 2: Increasing S(r) for anisotropic isothermal sphere; c = —0.5, k = 2 
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5 Example 2 



As a second example we demonstrate the applicability of the algorithm in by 
generating anisotropic Schwarzschild spheres. The line element for the interior 
Schwarzschild model [1] is 



ds 2 



(A - BA) 2 dt 2 + A~ 2 dr 2 + r 2 (d6 2 + sin 2 ddtf) 



(3? 



where A 
tions for 



1 — r 2 /R 2 } A and B are constants. The corresponding isotropic func- 
are given by 



(i/ , A ,m ,p ) 



2\n{A- BA}, -InA 5 



1 



2# 2 ' R 2 



A-3BA 
A-BA 



(39) 



The energy density function has the form /i = 3/R 2 . We therefore have 



/'o 



constant 



(40) 



for incompressible matter. This is a reasonable approximation in particular situ- 



ations as the interior o: 
of near uniform density 



dense neutron stars and superdense relativistic stars are 



24 



Consequently the assumption ()4()j) of uniform en- 



ergy density is often used to build prototypes of realistic stars in the modelling 



process[ 



The integral J a in (|2fij) takes the form 

2Br 2 



2Br 3 - r 1 



2B 2 r 



2^,2 



+ 



2B 



i? 4 A 3 (A + 5A) R*A 2 (A + BA) 2 R 2 A(A + BA) 



ABr 2 



R 2 A (A + BA) 



x 



i? 2 A (A + BA) (l + 
6 



2Br 2 



2 + 



2Br 2 



-i 



+ 



V R 2 A (A + BA) J \ R 2 A(A + BA)J_ 
2Br \ / 2Br 2 \ / 2Br 2 



r R 2 A (A - BA) 



1 + 



R 2 A (A - BA) 



2 + 



R 2 A (A - BA) 



(41) 



dr 
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With the substitution u = A = a/1 — r 2 /R 2 , (}4~Tj) becomes 

(B + 3Au - ABu 2 ) (2B + Au - 3Bu 2 ) 



J a 



[1 - u 2 ) (A - Bu) (B + Au - 2Bu 2 ) 
Wu 2 (A - Bu) 
~ (2B + Au- Wu 2 ) (B + Au- 2Bu 2 ) 

2B 2 (1 - u 2 ) (A - 2Bu - 2Au 2 + Wu 3 ) \ 
+ u (A - Bu) (2B + Au- Wu 2 ) {B + Au - 2Bu 2 ) ) U 

The above integral can be simplified with the help of partial fractions. We obtain 

'1 3 3 B A-6Bu 

J a = / - + r - 7TTT-. ; + ~. 7T~ + 



u 2{l-u) 2 (1 + u) A- Bu 2B + Au - Wu 2 
2A - 7Bu 



B + Au- 2Bu 2 



du 



3 3 

In u - - In {1 - u} - - In {1 + u} - In {A - Bu} + In {2B + Au - Wu 2 } 



A/B (7/2) (u - A) 



A 2 +8B 2 _ (.. _ _A_\ 2 A 2 +8B 2 _ ( _ A_\ 
WB 2 \ U 4Bl 16B 2 \ U AB) 

,2\ 



du 



In 



u (2B + Au- Wu 1 



3 



(A - Bu) (1 - u 2 Y (B + Au- 2Bu 2 Y 

We have completed the integration and obtained J a in terms of the intermediate 
variable u. In terms of the original variable r (in A) we can write J a as 

A (2B + AA- 3BA 2 ) 



J a = In 



(A - BA) (1 - A 2 f 2 (B + AA- 2BA 2 ) 7/A 



A t f 1 y/^+SB 2 ( A 4g) 

4V^+8^ n ll + 7 ^ 3f (A- 




Then the function a in (1261) becomes 



kR*A (2B + AA- 3BA 2 ) ( 1 - T^ftgf (A - j|) \ 4 vTO^ , 
r3 (A - BA) (5 + AA - 25A^ 4 1 + 73 |ft 5 , (A - A) _ I 
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and {2Zj) and 



x 



-InU 



respectively lead to 

kR 3 (2B + AA- 3BA 2 ) 
' r A (A - BA) {B + AA- 2BA 2 ) 7/4 V ' # 2 A ( A ~ BA ) 



2Br 2 



W (A -A) 



VA 2 +8B 2 



i\/A 2 +8B 2 



■IB 



VA 2 +8B 2 



(A -A) 



fci? 3 A (25 + AA - 3B A 2 



1 + 



25r 2 



2 (A — BA) (-B + AA - 2BA 2 ) 7/4 \ 5 2 A (A - BA) 



1 - 



IB 



VA 2 +8B- 



(A - j|) \ 4V^ 2 +8S 2 



1 + 



4B 



(A 



4Bj 



Hence the new line element has the form 



ds 2 



(A-BAfdt 1 



1 

+ A^ 



1 + 



kR 3 (2B + AA- 35 A 2 



1 + 



25r 2 



1 



x 



1 + 



rA (A - BA) (B + AA- 2BA 2 ) 7 ' 4 V ^ 2 A (A - BA) 

dr 2 + r 2 (d6 2 + sin 2 6d(f) 2 ) 



IB 



VA 2 +8B 



(A - fg) \ Wa 2 +sb 2 



'IB 



-1 



VA 2 +8B 2 

and the matter variables have the form 

r 3 LR 3 A (25 + AA- 3BA 2 ) 



m 



2R 2 



x 



2 (A - BA) (B + AA- 2BA 2 ) 7/4 \ 5 2 A (A - BA) 



1 + 



25r 2 



IB 



VA 2 +8B 2 



(A - f^) \ WA 2 +8B 2 



4B 



VA 2 +8B 2 



(A -A) 



A-3BA 



+ 



kR 3 A (2B + AA- 35 A 2 



R 2 (A - BA) r 3 _ # A ) (5 + AA - 25A 2 ) 7/4 



4B 
Vyl 2 +8B 2 

4£ 
VA2+8B 2 



(A - 4 V / I 2 >i^ 

W£) 
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Pi 



A-3BA 



kR 3 A (2B + AA- WA 2 



R 2 {A - BA) 2r 3 (A - BA) (B + AA - 2BA 2 ) 7/4 



1 



w 



VA 2 +8B- 



(A - ^) \ WA 2 +8B 2 



4,B 



(a -A) 



(4* 



1 VA 2 +8B 2 

The isotropic Schwarzschild sphere model ()38|) generates the anisotropic Schwarzschild 
sphere model ()45j1 - ()48j) . With the parameter value k = we regain the original in- 
terior Schwarzschild sphere. 



S 



The degree of anisotropy has the form 
\/3kR 3 A (25 + AA- WA 2 ) 



1 



■w 



VA 2 +8B 2 



(A-^)\*V^Sb 



(49) 



2r3 (A - BA) (B + AA- IB A 2 ) 1 " \\ + (A - A) 

Mathematical^ was once again used to plot the anisotropic factor ()49|) . The re- 
sulting plot is shown in Figure El for chosen particular values of the parameters A, 
B, k and R. Other choices of these parameters may produce a different behaviour 
for S. The plot of S against r is in the interval < r < 1. The fact that the 
anisotropic factor is in closed form and the profile as shown in Figure H3 shows that 
physical analysis of this model can be investigated, which will be pursued at a later 
stage. 



is 
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Figure 3: S(r) for anisotropic Schwarzschild sphere; A = 3, B 
R = l. 



1, k = 1 and 



14 



Acknowledgements 

SDM and MC thank the National Research Foundation of South Africa for financial 
support. MC is grateful to the University of KwaZulu-Natal for a scholarship. 

References 

[I] M S R Delgaty and K Lake, Comput. Phys. Commun. 115, 395 (1998) 

[2] M R Finch and J F E Skea, Preprint available on the web: 
:/ /edradour. symbcomp.uerj.br/pubs. html (1998) 

[3] H Stephani, D Kramer, M A H MacCullum, C Hoenslaers and E Her It, Exact 
solutions of Einstein's field equations (Cambridge University Press, Cambridge, 
2003) 

[4] K Dev and M Gleiser, Gen. Rel. Grav. 34, 1793 (2002) 

[5] K Dev and M Gleiser, Gen. Rel. Grav. 35, 1435 (2003) 

[6] L Herrera, J Martin and J Ospino, J. Math. Phys. 43, 4889 (2002) 

[7] L Herrera, A D Prisco, J Martin, J Ospino, N O Santos and O Troconis, Phys. 
Rev. D 69, 084026 (2004) 

[8] B V Ivanov, Phys. Rev. D 65, 10411 (2002) 

[9] M K Mak and T Harko, Chin. J. Astron. Astrophys. 2, 248 (2002) 
[10] M Mak and T Harko, Proc. Roy. Soc. Lond. A 459, 393 (2003) 

[II] R Sharma and S Mukherjee, Mod. Phys. Lett. A 17, 2535 (2002) 

15 



[12 
[13 
[14 
[15 
[16 

[it; 

[18 
[19 
[20 
[21 
[22 

[23 
[24 
[25 
[26 



S Rahman and M Visser, Class. Quantum Grav. 19, 935 (2002) 

K Lake Phys. Rev. D 67, 104015 (2003) 

D Martin and M Visser, Phys. Rev. D 69, 104028 (2004) 

P Boonserm, M Visser and S Weinfurtner, ArXiv:gr-qc/0503007 (2005) 



S D Maharaj and M Chaisi, Math. Meth. Appl. Set. 29 67 (2006) 

S D Maharaj and R Maartens, Gen. Rel. Grav. 21, 899 (1989) 

M K Gokhroo and A L Mehra, Gen. Rel. Grav. 26, 75 (1994) 

M Chaisi and S D Maharaj, Gen. Rel. Grav. 37 1177 (2005) 

M Chaisi and S D Maharaj, Pramana - J. Phys. Submitted (2006) 

W C Saslaw, S D Maharaj, and N Dadhich, Astrophys. J. 471, 571 (1996) 

W C Saslaw, Gravitational physics of stellar and galactic systems (Cambridge 
University Press, Cambridge, 2003) 

S Wolfram, Mathematica (Wolfram, Redwood City, 2003) 

S D Maharaj and P G L Leach, J. Math. Phys. 37, 430 (1996) 

C E Rhoades and R Ruffini, Phys. Rev. Lett. 32, 324 (1974) 

R L Bowers and E P T Liang, Astrophys. J. 188 657 (1974) 



16 



